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15 Abstract 

The natural variability in the terrestrial mesosphere needs to be known to correctly quantify 
global change. The response of the thermal structure to solar activity variations is an 
important factor. Some of the earlier studies highly overestimated the mesospheric solar 
20 response. Modeling of the mesospheric temperature response to solar activity has evolved in 
recent years, and measurement techniques as well as the amount of data have improved. 
Recent investigations revealed much smaller solar signatures and in some case no significant 
solar signal at all. However, not much effort has been made to synthesize the results available 
so far. This article presents an overview of the energy budget of the mesosphere and lower 
25 thermosphere (MLT) and an up-to-date status of solar response in temperature structure based 
on recently available observational data. An objective evaluation of the data sets is attempted 
and important factors of uncertainty are discussed. 
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INTRODUCTION 


The study of variation in atmospheric parameters due to several natural periodic and episodic 
events has always been an interesting subject. It was realized recently that the perturbation of 
35 atmospheric parameters caused by various human activities is not only confined to the lower 
atmosphere but that it also most likely extends into the upper atmosphere [Roble and 
Dickinson, 1989; Roble, 1995; Beig, 2000]. In view of this, it has become all the more 
important and vital to study the variations due to natural activities in parameters affecting 
climate to distinguish them from perturbations induced by global change. 

40 

Variations arising on decadal and even longer time scales may play a significant role in long- 
term trend estimates. One of the major sources of decadal variability in the atmosphere is the 
11 -year solar activity cycle (as modeled by Brasseur and Solomon, 1986). Electromagnetic 
radiation from the sun is not constant and varies mainly at shorter UV wavelengths on 
45 different time scales [Donnelly, 1991]. Incoming solar radiation provides the external forcing 
for the Earth-atmosphere system. While the total solar flux is quite constant, the UV spectral 
irrandiance on the timescale of the 27-day and 11-year solar cycles exhibits the largest 
changes, up to a factor of two over a solar cycle for the solar Lyman-a flux. Studies on the 
changes in solar UV spectral irradiance on timescales of the 27-day and 11 -year solar cycles 
50 have been attempted by many workers in the past [Donnelly, 1991, Woods and Rottman, 
1997, etc.]. It is believed that the essentially permanent changes arising in several 
mesospheric parameters due to human activities are weaker whereas periodic changes due to 
variations in solar activity are comparatively stronger [Beig, 2000]. 

55 The study related to the influence of solar activity on the vertical structure of temperature and 
its separation from global change signals has been a challenge because only data sets of short 
length (one or two decades) were available. The analysis of systematic changes in 
temperature in the mesosphere and lower thermosphere has not been as comprehensive as in 
the lower atmosphere. It is possible to suppress or even almost avoid the effects of solar cycle 
60 on trend determination with the use of proper selection of the analyzed period combined with 
the use of data corrected for solar and geomagnetic activity, or by comparison with empirical 
models, which includes solar and geomagnetic activity, local time, season, latitude and 
maybe some other parameters. The solar and geomagnetic activity may have a crucial impact 
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on the trend determination when data series are relatively short, or when we study trends in 
65 the ionized component (ionosphere). Modeling of the mesospheric data series to extract the 

solar cycle response has evolved with time, as improvements have been made in the 
measurement techniques of the 11 -year solar UV spectral changes. Most of the earlier 
predictions overestimated the mesospheric response, since they were based on incorrect solar 
UV-radiation derived from data of insufficient quality and/or length. This situation only 
70 changed with the SME and UARS-missions, when the data became available to quantify the 

variations since 1981. The modeling work of Chen et al. [1997] reported a solar cycle 

response of several Kelvin in the mesopause region. The observed temperature variability at 
70 km is not explicable in terms of corresponding 11 -year changes in observed ozone 
[Keating et al., 1987]. Searches for a strong dynamical feedback and attempts to invoke a 
75 strong odd hydrogen photochemical heating effect, have so far not been successful. Until 
recently, different data sets showed solar-cycle responses different even in polarity. The 
limited availability of data sets and the comparatively short length of data records have been 
the major constraints for mesospheric analysis. Nevertheless, during the past decade, a 
number of studies have been carried out and more reliable solar signals in mesospheric 
80 temperatures have been reported. It was thought earlier that it would be hard to identify a 
trend in the MLT region if solar response is very large in magnitude, and that we needed 
longer data sets encompassing several solar cycles. Recent investigations revealed the 
presence of a solar component in MLT temperature in several data sets but probably they are 
not as strong as expected. In recent times, a number of studies related to 1 1 -year periodicities 
85 in temperature of the MLT region have been reported [She and Krueger, 2004 and references 
therein]. 

However, the solar response in temperature, if not properly filtered out, is still one of the 
major sources of variation which may interfere with the detection of human-induced 
90 temperature trends for the MLT region [Beig, 2002] and will have strong implication in the 
quantification of global change signals. In recent time, the search for the effects of the 11 -yr 
solar cycle on middle atmosphere temperature has not led to consistent results that were easy 
to interpret. Model studies suggest an in-phase response to the UV flux, peaking in the upper 
mesosphere (2 K ampl.) and at the stratopause (1 to 2 K ampl.) [e.g., Brasseur, 1993; Matthes 
95 et al., 2004], However, the satellite analysis of Scaife et al. [2000] indicates a maximum 
response at low latitude of about 0.7 K between 2 and 5 hPa (around 40 km), while that of 
Hood [2004] shows a near zero response at 5 hPa but then increasing sharply to 2 K near 1 


3 



hPa. The increase of solar influence with altitude is not smooth. For example, the solar effect 
in the mesopause region is relatively small (according to the model by Matthes et al., 2004; 
100 but also according to several observations, see below), so it is easier to study long-term trends 
in this region. It should be clear that if one does not account properly for the solar cycle 
response, there can be biases for any remaining trend term. This concern is a particular 
problem for any data time series that is not well-calibrated, not representative of seasonal or 
global-scale processes, or not long enough. 

105 

Because of the very limited data not much effort has been made to synthesize the results 
available in the past. Consequently, our knowledge of the quantification of solar response in 
the temperature based on observations and model calculations for this region has been rather 
poor. In view of this, it would be highly desirable that a consolidated status report for solar 
110 trend in thermal structure for this region be prepared. 

Before ground-based instrumentation with sensitive photoelectric registration and rocket- 
borne in-situ measurements became available, the search for solar cycle effects with visual 
airglow photometry in the 1 920s, and with photographic spectrography, still dominant in the 
115 1960s, is now mainly of historical interest. Reliable temperature determinations, by whatever 

technique, became available only a few decades ago. Only recently, the detection of solar 
activity effects in the upper atmosphere comes close to become a routine affair, and the 
length of the available data sets is the main factor determining the quality of the results. In 
order to arrive at a balanced overview of our present knowledge, it is therefore natural to 
120 focus on the most recent results. These are also often based on the longest-duration data sets 
of homogeneous quality. The reader interested in the historical development can find 
references about early investigations of the atomic oxygen green line, which date back to the 
1920s, and of subsequent Doppler temperature determinations since the mid-1950s, in 
Hernandez [1976]. Other useful references focussing on green-line intensity variations can be 
1 25 found in Deutsch and Hernandez [2003]. 

This article reviews the present status of observational and modelling evidence on the 
response of the temperature structure in the region from 50 to 100 km to solar activity 
variations. An objective evaluation of the available data sets is briefly attempted and 
130 important uncertainly factors are outlined. We also discuss the lower thermosphere briefly. 
For convenience, the whole region from 50-100 km is referred to as the MLT (Mesosphere 
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and Lower Thermosphere) region. The region from 50-80 km will be referred to as the 
“mesosphere” and the region 80-100 km will be referred to as the “mesopause region”. 
Understanding and interpreting the causes of atmospheric trends requires a fundamental 
135 understanding of the energy budget. This is essentially the focus of the entire field of 
tropospheric climate science, which is seeking to determine the extent to which human 
activities are altering the planetary energy balance through the emission of greenhouse gases 
and pollutants. We are just now at the point of being able to quantitatively assess the energy 
budget in the MLT for the first time using the Thermosphere-Ionosphere-Mesosphere 
140 Energetics and Dynamics (TIMED) mission and Sounding of the Atmosphere using 
Broadband Emission Radiometry (SABER) instrument data. 

2. OVERVIEW OF THE ENERGY BUDGET OF THE 
MESOSPHERE AND LOWER THERMOSPHERE 

145 

Earth’s mesosphere and lower thermosphere are regions in which the transport and exchange 
of energy occur through subtle and complex processes. The main inputs to the system are of 
course provided by the Sun in the form of both photon and particulate energy. Ultraviolet 
(UV) radiation from 1 to 300 nm is absorbed primarily by molecular oxygen and ozone in the 
1 50 MLT. The variability of this portion of the spectrum with the 1 1-year solar cycle affects both 
the thermal structure (through changes in the overall amount of energy deposited) and the 
photochemistry of the MLT, especially the ozone abundance. Ozone is of particular 
importance to the MLT energy budget. Through the absorption of solar radiation and as a 
participant in exothermic chemical reactions ozone is responsible for up to 80% of the solar 
155 and chemical heating of the mesosphere [Mlynczak, 1997]. Here we will provide a brief 
overview of the energy budget of the MLT region in this section, following the corresponding 
presentation in Beig et al. [2003] where more details are given. 

The critical elements of the MLT energy budget are: heating due to the absorption of solar 
160 radiation by O 2 , O 3 , and CO 2 ; cooling due to infrared emission from NO, CO 2 , O 3 , H 2 O and 
O; heating due to exothermic chemical reactions involving odd-oxygen and odd-hydrogen 
species; energy loss due to airglow emission by 02 ('A), 02 ('E), 002 ( 4.3 pm) and OH(u). It is 
important to distinguish the energy loss due to airglow from that which is characterized as 
cooling. The energy in the airglow reduces the efficiency of solar or chemical reaction 
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heating - it never enters the thermal field, and hence is not acting to reduce the kinetic 
temperature of the atmospheric gases. Finally, particle input, especially in the thermosphere, 
is important, especially on short timescales and associated with solar flare or coronal mass 
ejection (CME) events. A summary of key heating (solar and chemical) and infrared cooling 
terms is given in Table 1. 

170 

The single most significant factor in differentiating the energy balance of the MLT from the 
atmosphere below is that the density in the MLT is so low that collisions cannot always 
maintain the processes of absorption and emission of radiation under local thermodynamic 
equilibrium (LTE). Consequently, computation of rates of solar heating and infrared cooling 
175 is a much more challenging process. In the case of solar energy, not all of the absorbed 
energy is thermalized locally. This fact requires a detailed accounting of all possible 
pathways for the absorbed solar energy to transit prior to ending up as heat or to being 
radiated from the atmosphere as airglow without ever having entered the thermal field. Thus 
we say that the efficiency of solar heating is substantially less than unity due to these 
1 80 processes, to as low as 65%. The details of the solar and chemical heating and the associated 
efficiencies are reviewed by Mlynczak and Solomon [1993]. 

In the case of radiative cooling, the effective temperature at which infrared active species 
radiate is not given by the local kinetic temperature. This fact requires extremely detailed 
1 85 consideration of the exchange of energy (thermal, radiative, chemical) between the infrared 

active molecules and their environment for a multitude of quantum energy states within each 
molecule [e.g., Lopez-Puertas and Taylor, 2002], The key radiative cooling mechanisms in 
the MLT involve several infrared active species including the molecules CO2 and O3 [Curtis 
and Goody, 1956], H2O [e.g., Mlynczak et al., 1999], NO [Kockarts, 1980], and atomic 
190 oxygen (O) [Bates, 1951]. Of these, the CO2, O3, and NO emissions exhibit substantial 

departure from LTE in the MLT. The water vapor and atomic oxygen emissions correspond 
to transitions in the far-infrared portion of the spectrum (wavelengths typically longer than 20 
pm) that are more readily thermalized by collisions, and thus maintained in LTE. 

195 The solar photon energy is the dominant source of energy into the MLT, but the solar 

particulate energy is nevertheless important. While the photon energy from the Sun varies on 
relatively long timescales (from the 27-day solar rotation to the 11 -year sunspot cycle), 


6 



particulate energy from the Sun varies in a much more erratic (and often violent) way. Recent 
observations of the thermospheric and mesospheric response to variations in particle input 
200 from CME events clearly indicate the potential to alter the thermal structure and the radiative 
cooling mechanisms. Seppdld et al., [2004] and Rohen et al., [2005] have observed the 
destruction of ozone in response to strong solar storm events. In the stratosphere and lower 
mesosphere, radiative cooling by ozone is critical to the energy balance. Thus there is a direct 
impact on the energy balance in the stratosphere from solar particle precipitation. 

205 

From this overview of the energy budget, it is clear that the variability of solar radiation input 
into the MLT region may impact the thermal structure directly (through the increase or 
decrease in the total amount of solar energy deposited) and indirectly, for example, by 
modifying the ozone abundance and thereby the heating and cooling rates. It is specifically 
210 because of the complexity of the energy budget that assessing and attributing observed 
changes (cyclical and secular) in the MLT is a formidable scientific challenge. 

3. TECHNIQUES AND OBJECTIVE EVALUATION OF 
TEMPERATURE DATA SETS 

215 

In addition to the satellite datasets, there are a number of experimental data records from 
ground based or in-situ observations of mesospheric temperatures, although the mesospheric 
record is still small compared to what is available at lower altitudes. The same techniques and 
associated measuring uncertainties that are discussed by Beig et al. [2003] for their use in 
220 trend analyses are also relevant, to some extent, in the present context, so we can make 
reference to the greater detail given there. 

In addition to ground-based observations which are capable to supply long data sets at fixed 
geographic locations, also in-situ data from rocketsondes and global observations from 
225 satellites can be used to measure temperature suitable for the detection of solar activity 
effects. Details of all the data sets obtained during the past few decades and available for 
evaluation of temperature trend in the MLT region are also given in Beig et al. [2003]. Figure 
1 shows most of the known ground based locations of long term temperature measurements 
all over the globe. The techniques applied to measure the temperature are also indicated as far 
230 as possible. As mentioned in Beig et al. [2003], even for standard instrumentation used for a 
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long time, technical improvements can introduce uncertainties when data obtained at different 
times are combined into longer data sets. While this can be most serious for long-term trend 
detection, it can also interfere to some extent with the determination of solar activity 
response. Techniques capable of supplying data over an extended height range like rocket- 
235 launched or lidar temperature or density soundings (from which temperature profiles can be 
derived) nevertheless suffer from an inevitable loss of precision at the greatest altitudes, 
where they often cannot compete with ground-based observations. The long time span 
covered by some ground-based measurements makes them particularly useful to study solar 
cycle variations. 

240 

Mesopause region temperature is most often determined from line intensity measurements in 
hydroxyl (OH) airglow bands of the airglow, but the so-called Atmospheric band of 
molecular oxygen is now quite often also used. These rotational temperatures agree with 
kinetic temperature at the peak of the vertical airglow emission profile. According to 
245 measurements with many different techniques, the OH emission comes from an emission 
layer at 87 km with a mean thickness of 8 km [Baker and Stair, 1988; more references are 
given in Beig et al., 2003]. Satellite limb scans have resulted in reports on height variations 
by several kilometers which may be related to dynamics [e.g., Liu and Shepherd, 2006], and 
Nikoukar et al. [2007] have found recently that the bands from the upper vibrational levels 7, 
250 8, and 9 come from an altitude slightly higher than the bands from the 4, 5, and 6 levels. 

According to these results, the difference is 1 .9 ± 1 .4 km, while the mean peak altitude for the 
latter bands (which are probably representative of the most widely used ground-based 
observations) is consistent with the nominal values mentioned above. The observed 
variability does not invalidate ground-based measurements of hydroxyl rotational 
255 temperature as a useful tool to diagnose atmospheric temperature trends or solar activity 
effects, as long as this variability can be treated as random, or be considered as part of the 
phenomenon. The same holds for the O 2 Atmospheric band, with a nominal emission peak 
height of 95 km. 

260 From OH or O 2 airglow observations, temperature precisions of a few Kelvin can be obtained 
with integration times not longer than a few minutes. Therefore, by averaging over a number 
of individual measurements, the contribution of instrumental noise to the mean temperature 
can easily be made negligible. Systematic errors affecting data accuracy have only an 
influence on trend or solar activity results if they vary with time. They are not a problem if 
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265 long-term stability can be assured, and one way to ensure stability is with good instrument 
calibration. The discussion about this point in Beig et al. [2003] is mostly important in the 
context of the possibility to detect small long-term trends. For detecting effects of the 1 1 -year 
solar cycle, which has a rise time of only about four years, and where responses of several 
Kelvin have been reported, the instrument stability requirements are less stringent. The 
270 relative calibration of the instrument response at two or more wavelengths necessary for 
determining rotational temperature is not difficult. 

The rotational temperatures in the mesopause region vary on time scales from a few minutes 
for short-period gravity waves to the solar cycle, and beyond. Nocturnal mean temperatures 
275 used as the basis for solar cycle and trend analysis are affected by the short-term variability 
only as far as gravity waves and variations due to the thermal tide are not completely 
cancelled out. This “geophysical noise” can be expected to be quite variable and so create 
only small uncertainties on longer time scales. There is however also a day-to-day variability 
from planetary waves and unknown sources which could not be avoided even if 
280 measurements over complete nights were always available. This underlines the importance of 
dealing with airglow temperature data sets based on the greatest possible number of nights. 
The same obviously also holds for data sets from other techniques. 

Apart from O 2 rotational temperatures, some data sets extend the information available above 
285 the altitude of OH by using atomic line intensities from sodium or atomic oxygen as a proxy 
for temperature, based on an empirical correlation between intensity and temperature (see, 
e.g., Golitsyn et al., [2006]). The validity of this approach is questionable and cannot be 
recommended as a replacement for direct temperature measurements, be it by the 
measurement of O 2 rotational temperature, of Doppler width with Fabry-Perot instruments, or 
290 by laser spectroscopy with sodium lidars. 

The length of the data set required for determining solar signatures may be as short as the few 
years that the solar cycle takes to ascend from minimum to maximum, but may also be as 
long as several cycles, if the effect is small compared to other variability (for example, 
295 seasonal) , or is itself strongly variable. The data sets from different measurement techniques 
vary widely not only in the number of years covered, but even more so with respect to the 
uniformity of coverage and the number of individual data points available. Some airglow data 
sets consist of millions of individual, statistically independent observations at a fixed site, 
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resulting in up to about 5000 nocturnal means, all referring to the same (nominal) altitude. On 
300 the other hand, rocket soundings yield only one profile per launch, and from less than 100 to 
several hundred profiles may be available from a given site, but a considerable altitude range 
is covered. As pointed out before, lidar soundings (either by Rayleigh lidars covering a wide 
range of altitudes similar to rocketsondes, or by sodium lidars that are limited to the 
mesopause region) can easily surpass the number of profiles from rockets, being limited only 
305 by clear weather requirements, and not by equipment expense. Finally, satellite observations 
easily comprise millions of vertical temperature profiles, with near-global coverage, but the 
number of overpasses at a given site is very much lower, refers to only slowly varying local 
time, and the available long-term coverage is still small. The SCIAMACHY instrument on 
Envisat is capable of measuring OH rotational temperature by limb sounding [von Savigny et 
310 al., 2004], It was launched only in 2004, but can be expected to contribute data on solar 
activity response in the near future. 

4. OBSERVATIONS OF THE SOLAR RESPONSE IN THE 
315 MESOSPHERE 

The response of temperature and the middle atmosphere species to 11 -year solar UV 
variations has been difficult to isolate using satellite data. This is partially due to the short 
time series of satellite data sets relative to 1 1-year variations, to instrument drifts, and to the 
320 strong longitudinal variability that makes zonal means appear quite noisy [see Chanin, 2006]. 
On this time scale the quasi coincidence of the recent major volcanic eruptions with solar 
maximum [Kerzenmacher et al., 2006] conditions increases the challenge while indirect 
mesospheric responses were observed [Keckhut et a.l, 1995, 1996]. This effect was caused by 
changes to the wave propagation induced by the thermal forcing inside the volcanic cloud and 
325 vertical stability around the tropopause [Rind et al., 1992], In the past, only rocket 
temperatures have provided such long data sets in the mesosphere. However, the required 
aerothermic corrections and changes of the sampling and time of measurements induce some 
bias mainly in the mesosphere. More recently Rayleigh lidars that are much less expensive 
and require less resources for continuous operations have replaced the rocket techniques. 
330 From space, HALOE aboard UARS is the only experiment that allows mesospheric 
temperature over more than a decade with a single instrument. However, while global, the 
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number of solar occupation does not provide such a large sampling as desirable. 


The solar activity is also modulated by the solar rotation, and UV series exhibit strong 
335 responses with periods of 27 days and harmonics. On this scale more temperature series are 
available. From the ground, rockets and lidars can be used. However, lidars are more 
adequate to perform daily series while typical rocket sampling is close to a week. The 
experiments on board Nimbus 6 and 7 were used intensively to retrieve stratospheric and 
mesospheric temperature responses. However, these changes related to the solar rotation 
340 present smaller amplitudes than the solar cycle and are highly non-stationary. In the case of 
the sun, the physical processes governing the evolution of active regions and the resulted 
variations in the solar output are, at best, only quasi-stationary over a limited time period. 

On the time range of solar cycle the radio flux at 10.7 cm is used as proxy of solar activity 
345 while long-term UV measuremeots from space are not available. On the other hand, the 
short-term solar UV variation is not well described by standard radio solar flux at 10.7 cm 
and more direct UV measurements from space at 205 nm, Lyman alpha or proxies such as 
Magnesium Lines Mg II are preferred to better describe daily changes of solar UV (see 
Dudok de Wit et al., [2007] for a recent investigation of this topic). While in the middle 
350 atmosphere, ozone and temperature are highly connected due to thermal ozone absorption 
and thermal sensitivity of the ozone dissociation, the simultaneous investigations of ozone 
and temperature allow for a better understanding of the middle-atmosphere response to solar 
activity changes. Ozone measurements on Solar Mesosphere Explorer (SME) and on Nimbus 
7 were analyzed. At high latitudes, direct ozone response (and hence temperature) to solar 
355 activity variations may also be overwhelmed by solar proton events. The satellite sensors for 
solar EUV may also occasionally be saturated by solar particles [e.g. see an example in 
Scheer and Reisin, 2007]. 

In photochemical models, the ozone sensitivity on the 27-day and the 11 -year scale are 
360 similar because the time constant of ozone is negligibly short, in comparison. However, 
discrepancies exist when including temperature-chemistry feedback in the model 
calculations. It is possible that this is indicative of an indirect dynamical component of the 
solar response. 
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4.1 Changes due to 27-Day UV Solar Forcing 


Temperature variations are affected by a number of short-term dynamical influences 
independent of solar variations, and thus it is more difficult to isolate the solar signal. 
Temperatures are available simultaneously from the SAMS instrument on the Solar 
370 Mesosphere Explorer (SME) satellite. In the ±20° latitude band at 2 mbar a temperature 
variation of 1.5K for 10% ozone change is reported, which grows to 2.5K at 70 km. In 
contrast to the stratospheric maximum that is limited to the ±20° latitude band, this second 
maximum in the mesosphere is present in the ±40° latitude band. The observed temperature 
phase lag with 205 nm solar flux is shorter than 1 day in the mesosphere, and the altitude of 
375 maximum temperature sensitivity is close to the altitudes of maximum ozone depletion. 
Therefore, in addition to the HO x effect on ozone, the temperature sensitivity can be expected 
to play a role through temperature feedback or as a consequence of the solar Lyman alpha 
heating. This mesospheric maximum was not predicted by numerical models but must be 
real, since Summers et al. [1990] conclude that the discrepancies between model and 
380 observation cannot be explained by data-related errors. At mid-latitudes, wavelet analysis of 
lidar time series [Keckhut and Chanin, 1 992] shows that planetary waves tend to mask the 
direct solar response in temperatures since wave amplitudes are large and periods may be 
comparable to the solar rotation period; planetary waves exhibit periods from a few days to 
30 days. On the other hand, planetary waves may be directly involved in the solar forcing 
385 (see below). 

Ebel et al. [1986] have suggested that the generation of vertical wind oscillations in the 27- 
day period range would at least lead to the right sign of correlation (through adiabatic cooling 
during the upward wind phase and simultaneous transport in the direction of the vertical 
390 ozone mixing ratio gradient in the lower dynamical regime and photochemical increase at 
higher layers due to the temperature dependence of the ozone reaction coefficients). This 
effect may also be responsible for the fact that the ozone perturbations inferred from UV flux 
changes are better reproduced by simulations without temperature feedback than with it 
[Keating et al., 1987] due to the compensating effect of adiabatic heating. 

395 

Radiance measurements were made with the Pressure Modulated Radiometer (PMR) on 
board NIMBUS 6 [Crane, 1979]. Maximum values obtained for the 27 day periodicity were 
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1 .5 K near the mesopause, 3.0 K in the lower mesosphere and 3.5 K in the upper stratosphere, 
at latitudes between about 50° and 70° [Ebel et al., 1986]. Since indirect perturbations seem 
400 to exceed the direct ones in amplitude, non-linear interactions of forced variations with the 
atmospheric system also have to be considered. Furthermore the large spatial scales of the 
possible solar activity effects showing up in the global and hemispheric data sets employed in 
this study support the view that planetary waves are an essential part of the unknown 
mechanisms. 

405 

4.2 Changes on the 11 -Years Solar Scale 

The atmospheric temperature response to solar UV changes is expected through ozone and 
oxygen absorption processes in the middle atmosphere. In the stratosphere, the response 
410 shows a positive correlation between the temperature and the solar activity with an effect of 
1-2 Kelvin in the upper stratosphere due to ozone photolysis and solar absorption, while at 
higher latitude negative responses are reported [Keckhut et al., 2005]. These observations 
confirmed by rocket series could be explained by dynamic feedbacks and more specifically 
by the occurrence of stratospheric warmings [ Hampson et al., 2005]. From this numerical 
415 simulation, a positive effect is expected in the mesosphere. Because the winter response 
results in a dynamic feedback, the signature is expected to be non-zonal in the northern 
hemisphere [Hampson et al., 2006]. While stratospheric warmings are associated with 
mesospheric cooling, it is not surprising to see these alternating patterns at mid and high 
latitudes [Matsuno, 1971]. In the tropical mesosphere, a response can be also expected, as 
420 tropical mesospheric anomalies associated with stratospheric warmings are also reported 
[Sivakumar et al., 2004, Shepherd et al., 2007]. A summary of the solar response in the 
mesosphere is given in Table 2. 

The search for a solar trend in the mesosphere had started in the late 70s when a few authors 
425 reported solar cycle associated variability in mesosphere temperatures. Shefov [1969] 
reported a solar cycle variation in OH rotational temperature on the order of 20-25 K for mid- 
latitudes. Labitzke and Chanin [1988] using rocketsonde data at Heiss island located at 81°N 
reported solar cycle temperature variations on the order of 25K at 80 km. Kubicki et al. 
[2007] have reanalyzed the same set of data and deduced a negative solar response of several 
430 kelvins. The time series of Russian rocketsonde measurements at four different sites 
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(covering low to high latitudes) revealed a substantial positive solar response in the 
mesosphere [Mohanakumar, 1985, 1995]. However, in recent time, results are found to be 
quite different. The reanalysis of the Thumba (8°N) tropical data extending to more than two 
solar cycles [Beig and Fadnavis, 2003] has recently also resulted in a positive solar response 
435 of temperature in the mesosphere but of much lower magnitude than reported earlier [Kokin 
et al., 1990]. US rockets have shown a clear solar response in the upper stratosphere 
[Dunkerton et al., 1998]. In the mesosphere only 2 subtropical sites allow to retrieve the solar 
response in the mesosphere. A positive correlation has been found with a temperature 
response of 2 K on a large latitude range from 50 to 70 km [Keckhut et al., 1999].. An 
440 analysis of falling sphere and rocket grenade data by Liibken [2000, 2001] revealed no 
statistically significant solar component for the altitude range 50-85 km, however the analysis 
only included data during the summer season. Remsberg and Deaver [2005] have analyzed 
long-term changes in temperature versus pressure given by the long time series of zonal 
average temperature from the Halogen Occultation Experiment HALOE on the Upper 
445 Atmosphere Research Satellite (UARS). The HALOE temperature data are being obtained 
using atmospheric transmission measurement from its C02 channel centered at 2.8 
micrometers [Russell et al., 1993, Remsberg et al., 2002], While the length of the data set is 
still short, they have reported a mesospheric response of 2-3 K around 70-75 km. In a more 
recent work Remsberg [2007] found more accurate results ranging from 0.7 to 1 .6 K in the 
450 lower mesosphere and from 1.7 to 3.5 K in the upper mesosphere. In mid-latitude responses 
are larger and in the tropical latitude band only 0.4 to 1.1 K is reported. The long-term series 
of lidar data obtained at Haute Provence (44°N) has revealed a positive (in phase) solar 
response of 2 K/100 sfu (solar flux units) in the mesosphere up to 70 km. The response was 
found to fall off with height above 65 km, with a tendency towards a negative response above 
455 80 km [Keckhut et al., 1995]. 

Atmospheric temperature response to the 27-days and 11 -years solar cycle as function of 
altitude from different studies focus on tropical regions are shown in Figure 2. SAMS 
[Keating et al., 1987] and PMR [Ebel et al., 1986] analyses are performed on the solar 
460 rotation time scale. HALOE analyses above the ±20° latitude band [Remsberg, 2007] concern 
the 1 1-years time scale. Two different analyses [Beig and Fadvanis, 2003; Kokin et al., 1990] 
have been used to analyze the data of the rocket station Thumba (8.5°N, 77°E).. 
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Atmospheric temperature response to the 1 1 -years solar cycle as a function of altitude from 
465 different studies covering mid-latitudes are shown in Figure 3. Error bars associated with the 
OHP lidar correspond to the range of seasonal changes [Keckhut et al., 2005]. The hatched 
area corresponds to the HALOE observations in different latitude bands [Remsberg, 2007]. 
The US Rocket corresponds to an average of 3 subtropical rocket stations [Keckhut et al., 
1999]. The Soviet rocket sites (Heiss, Volgograd and Molodezhnaya) data have been 
470 analyzed by Kokin et al. [1990] and Golitsyn et al. (2006). They obtained different solar 
response for different stations as shown in Figure 3. 

4.3 Seasonal Variations 

475 Mohanakumar [1995] shows the summer response varied in the same way with latitude 
between Arctic and Antarctic but was about half the wintertime values. The Thumba results 
indicate also a stronger positive solar component in winter as compared to summer for the 
mesosphere, which is in agreement with mid-latitude lidar results [Keckhut et al., 1995]. 
Hauchecorne et al. [1991] had already reported earlier that solar response changes sign from 
480 winter to summer depending on height, using lidar data (1978-1989) at heights from 33 up to 
75 km. They found about 5 K/lOOsfu during winter for 60-70 km altitude (where the 
maximum response is observed) and for summer about 3K/100sfu. Later, Keckhut et al. 
[1995] using data from the same lidar for an extended period reported the solar response over 
a height range of 30-80 km for summer, with a negative tendency at height above 75 km. In 
485 the mesopause region, the changes of the response to solar activity during the specified 
intervals occur most distinctly in autumn, winter and spring. For summer, the response of the 
atmosphere does not practically change, though it has the maximal value. Therefore the 
dispersion of the values mentioned above is probably caused by the seasonal nature of 
observation. Golitsyn et al. [2006] have recently analyzed the response of the monthly mean 
490 temperature data on the solar activity variations for the altitudes 30-100 km. They obtained 
the minimal solar response at heights -55-70 km with a value of +2±0.4 K/lOOsfu for winter 
and -1±0.4K/100 sfu for summer. 

495 
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4.4 Atmospheric Response due to Solar Flares 


During major solar flare events, energetic particles penetrate down into the Earth’s 
mesosphere and upper stratosphere. By ionizing molecules, solar proton events (SPE) are 
500 expected to produce a large enhancement of odd nitrogen at high latitudes in the mesosphere 
[Crutzen, 1975; Callis et al., 2002], Odd nitrogen species play an important role in the 
stratospheric ozone balance through catalytic ozone destruction. In the upper stratosphere and 
mesosphere, ozone decreases of 20 to 40% associated with SPE have been reported [ Weeks et 
al., 1972; Heath et al., 1977; McPeters et al., 1981; Thomas et al., 1983, McPeters and 
505 Jackman, 1985]. When a strong stable polar vortex forms, diabatic descent inside the vortex 
can transport NO x rapidly downward, and may enhance the effect of ozone destruction 
[Hauchecorne et al., 2005]. As expected from numerical models [Reagan et al., 1981; Reid 
and Solomon, 1991] simultaneous cooling of around 10 K in the lower mesosphere was 
observed by Zadorozhny et al. [1994], in October 1989, using meteorological M-100B 
510 rockets, while NCEP temperature analyses [Jackman and Me Peters, 1985] report no 
detectable temperature decrease associated with the event of July 1982. 

A recent search for a response of the mesopause region to solar flares and geomagnetic 
storms by Scheer and Reisin [2007] in the airglow data base from El Leoncito (32°S, 69°W) 
515 revealed no convincing evidence, in spite of the coverage of very strong solar and 
geomagnetic events. There is however the possibility that the atmospheric response (at the 
relatively low latitude) is only short-lived and therefore limits detection by nocturnal 
observations to cases of favorable flare timing. If this is so, daytime observation techniques 
would be more suitable for detecting flare effects in the mesopause region. 

520 

5. OBSERVATIONS OF THE SOLAR CYCLE RESPONSE IN THE 
MESOPAUSE REGION 

5.1 Annual mean response 

525 

The majority of the results discussed in this section are obtained from measurements of OH 
airglow which correspond to nominal altitudes of 87 km. Section 4 of the Beig et al. [2003] 
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530 

535 

540 

545 

550 

555 

560 


paper contains information on some earlier solar activity studies since the 1970s, with an 
emphasis on OH rotational temperature, only part of which will therefore again be discussed 
here. 

Sahai et al. [1996] reported a solar activity effect of 32 K/lOOsfu from OH airglow 
temperature measurements at Calgary (51°N, 114°W) based on the comparison of two low 
activity years (1987, 1988) with a high activity year (1990). The solar signal found by 
Gavrilyeva and Ammosov [2002] from their OH data obtained at Maimaga (63°N, 130°E) 
between 1997 and early 2000 was only one third as great as the result from Calgary, but is 
still a high-end result (11+3 K/lOOsfu), when compared to the other solar temperature 
responses observed all over the globe, that were published in recent years. Both results may 
be considered problematic because of the short time span covered, in combination with the 
strong seasonal variability, characteristic of medium and high latitudes. They are not 
automatically discredited by disagreeing with the lower values found elsewhere, because they 
could only be refuted by observations under the same conditions. Similar arguments can be 
used in favor of the strongly negative solar cycle effects in O I Doppler temperatures (ca. -30 
K/lOOsfu) by Hernandez [1976] at mid- and low latitudes, and by Nikolashkin et al. [2001], 
at Maimaga (ca. -15 K/lOOsfu), although those latter authors’ result for OH temperature (+5 
K/lOOsfu for eastward QBO phase, -30 K/lOOsfu for westward) is somewhat at odds with the 
more recent data by Gavrilyeva and Ammosov [2002] from the same site. 

Figure 4 and 5 show the solar response in temperature (K/lOOsfu) as reported by different 
authors using various experimental data in recent time for the Northern and Southern 
Hemispheres, respectively. For easier reference, the pertinent details corresponding to each 
result are also listed in Table 3. 

Sigernes et al. [2003] found no solar signal in the time series of OH airglow data from the 
auroral station Adventsdalen (78°N, 15°E) that span 22 years. On the other hand, in the mid- 
latitude Northern Hemisphere, where the greatest number of OH airglow temperature 
measurements is available, all studies signal a positive response to solar activity. Espy and 
Stegman [2002] have initially not reported an appreciable solar cycle effect at the height of 
the OH layer over Stockholm (59°N, 1 8°E), but after a new analysis that includes more recent 
data, these authors find evidence for a positive effect of about 2 K/lOOsfu (1.6 + 0.8 
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K/lOOsfu, in winter); [P.J. Espy, 2007, personal communication]. New results from 
Zvenigorod (56°N, 37° E), for the years 2000 - 2006 reported by Pertsev and Perminov 
[2007] indicate an annual mean response of 4.5 + 0.5 K/lOOsfu, which is somewhat stronger 
than previous results partly derived from the same site [Golitsyn et al., 2006, see below]. 

565 The recent analysis by Offermann et al. [2004] based on 21 years of OH airglow temperature 
data for Wuppertal (51°N, 7°E) extends up to 2002 and now covers almost two solar cycles. 
This long series of observations was started in 1980 with the aim to determine solar and long- 
term trends in the mesospause region. The authors found an effect of 3.0 + 1.6 K/lOOsfu, on a 
monthly basis, from temperature enhancements during the maxima of two solar cycles. The 
570 authors assumed a linear correlation between temperature and solar activity, ignoring 
possible lags. The annual mean response is 3.4 K/lOOsfu. 

Bittner et al. [2000] analyzed the Wuppertal (51°N, 7°E) data for the period 1981-1995 with 
respect to temperature variability with periods of several days, but not with respect to 
575 absolute temperature (as most other studies). They found positive solar correlation response 
for temperature oscillations with periods greater than ~30 days and negative correlation for 
periods less than ~10 days. In an analysis of the complete Wuppertal data set including the 
year 2005 [Hoppner and Bittner, 2007] , the 1 1-year solar signal for the period range from 3 
to 20 days had disappeared, but there was, surprisingly, evidence for a correlation with the 
580 22-year heliomagnetic cycle. The authors investigated the possibility of a magnetic coupling 
between the solar and terrestrial magnetic fields, and indeed found a weak modulation in the 
Earth’s rotation period with a shape similar to the observed variation of the standard 
deviation of OH temperature. Lowe [2002] has made OH layer measurements for one full 
solar cycle at Delaware Observatory (43°N, 81°W) in Canada and found a positive solar 
585 response of 1.5 ± 1.1 K/lOOsfu. 

Golitsyn et al. [2006] have consolidated Russian results from their earlier analysis [ Golitsyn 
et al., 1996] based on the rocket data from Volgograd (49°N) and data from different airglow 
emissions obtained at Abastumani (42°N) and Zvenigorod (56°N), covering two activity 
590 cycles (1976-1991). Monthly and annual mean model profiles of the temperature response to 
the solar cycle where fitted to these data, which at the altitude of 87 km are based essentially 
on OH rotational temperatures. The annual mean response is about 2.7 +1.7 K/lOOsfu, with a 
tendency to grow with altitude. 
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595 An alternating negative and positive temperature response is consistently found in Northern 
Hemisphere mid-latitude results obtained by incoherent scatter radar, Rayleigh lidar [Chanin 
et al., 1989], and sodium lidar [She and Krueger, 2004] for altitudes between 30 and (in some 
data sets as high as) 140 km. This suggests that dynamical coupling from the troposphere to 
the thermosphere is involved in solar activity induced signatures. She and Krueger [2004] 
600 have recently reported the impact of 11 year solar variabilities on the mesopause region 
temperature over Fort Collins (41°N, 105°W) using their sodium lidar data obtained between 
1990 and 2001. They found no solar signal at 83 km, but a positive effect of 3 + 2 K/lOOsfu 
at 87 km (updated numbers according to C.-Y. She, personal communication [2007]). The 
HALOE results by Remsberg and Deaver [2005] mentioned above also refer to the lower 
605 limit of the mesopause region. For different latitude zones, 1 1-yr solar cycle terms with 
amplitudes of 0.5 to 1 .7 K were found for the middle to upper mesosphere (80 km). 

Figure 5 shows the solar response in temperature for the Southern Hemisphere. Clemesha et 
al. [2005] reported the OH rotational temperature measurements made at Cachoeira Paulista 
610 (23°S, 45°W) for the period from 1987 to 2000. A simultaneous linear and 1 1-year sinusoidal 

fit resulted in a solar cycle amplitude of 6.0 + 1 .3 K with maxima in 1990 (and therefore also 
in 2001), well in phase with solar activity. The linear trend of 10.8 + 1.5 K / decade agrees 
perfectly well with the results obtained for nearly the same time span and geographic area (El 
Leoncito, 32°S, 69 °W) by Reisin and Scheer [2002], Clemesha et al. [2005] also found a 
615 positive OH intensity trend of the order of 1.9%/year, which, in view of the combined error 
bounds, is not considerably above the intensity trend observed by Reisin and Scheer (about 
+l%/year). However, the strong solar cycle signature found by Clemesha et al. [2005] 
(expressed as 1 1-year amplitude) that can be estimated to correspond to about 8-10 K/lOOsfu 
is at odds with the near-zero effect enountered by Reisin and Scheer [2002], The Leoncito 
620 results are confirmed by the most recent analysis of OH (6-2) rotational temperatures and 
airglow brightness variations during the rise and maximum phase of solar cycle 23 [Scheer et 
al., 2005]. No solar cycle signature was found, but when a temporal trend is allowed for, the 
solar effect may approach 1.4K. The disagreement with the solar response at Cachoeira 
Paulista may be a consequence of latitudinal differences in planetary wave activity, and 
625 therefore need not be considered contradictory. 
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Results about OH temperatures from Davis (69°S, 78°E) will be discussed in the next section. 
Hernandez [2003] measured the polar mesospheric temperature above South Pole (90°S) 
from 1991-2002 using the Doppler width of the OH line at 840 nm by means of a high- 

630 reolution Fabry-Perot interferometer and deduced a solar signal as high as 13K/100sfu. This 

is the strongest solar temperature signal reported in recent years. From the same site, Azeem 
et al. [2007] have reported OH rotational temperatures obtained during the austral winters of 
1994—2004. In spite of the temporal overlap between both data sets, the comparable coverage 
of data, and the expected approximate equivalence of the temperatures obtained by both 

635 techniques, the solar cycle effect of 4.0 ±1.0 K/sfu reported by Azeem et al. [2007] is only 

about one third of the result obtained by Hernandez [2003]. The reason for this discrepancy is 
unknown. Figure 6 provides the summary of the solar response in temperature for the lower 
thermosphere. In this height range, the results obtained by Golitsyn et al. [2006] depend in 
part on temperatures estimated from the brightness of different airglow emissions, and in part 
640 on the extrapolation of OH temperature response to greater heights. Their composite profile 
indicates an annual mean solar response of 2.8K/100sfu and 2.3K /1 00 sfu for the altitudes 
of 92 km and 97 km respectively. She and Krueger [2004] have found a solar signal of 
4K/100sfu between about 92 and 98 km. Above 101 km, the effect decreases quickly to zero 
at 103 km and becomes negative at 104 km. The mean solar cycle effect in O 2 rotational 
645 temperatures measured at El Leoncito (32°S, 69°W) [Scheer et al., 2005] is consistent with 
the range of upper limits estimated earlier [Reis in and Scheer, 2002], A positive response of 
3.3 ± 0.3 K/lOOsfu is reported (which would reduce to 1.32 ± 0.3 K/lOOsfu, if a temporal 
trend is simultaneously fitted). The authors conclude however that the mean values are 
however only the net effect of successive short-term spells of anti-correlation, the absence of 
650 correlation, each lasting many months, and a 32-months regime of strong correlation. So, 
there is obviously no seasonal regularity in the solar signal at this site. All these recent results 
for northern and southern mid-latitudes that refer to heights around 95 km compare quite well 
with each other, although they are obtained by quite different techniques. 

655 5.2 Seasonal Differences 

As mentioned above, the results reported by Golitsyn et al. [2006] are seasonally resolved, 
and solar response changes with season (even in sign) and assumes the most extreme values 
(both positive, and negative) in the mesopause region, at heights of about 80-95 km. 
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660 Response is most variable in autumn, winter and spring, and a strong, but stable response 
prevails in summer. These authors deduced a solar response of -5+1.7 K/lOOsfu for winter, 
and 8+1.7 K/lOOsfu for summer, in this altitude range. 

The absence of a solar response in the OH data from Adventsdalen (78°N) for winter was 
665 already mentioned [Sigernes et al., 2003, and also Nielsen et al., 2002], Lubken [2000] 
arrived at a similar conclusion about 80 km at Andoya (69° N), from the comparison of 
rocket soundings, but those were made mostly during summer. Since most of the soundings 
correspond to low solar activity, this evidence is however not very strong. For high latitudes 
in the Southern Hemisphere (Davis, 69°S, 78°W), French and Burns [2004] have reported a 
670 positive solar response of 6 to 7K/100sfu in mid-winter, but smaller values (even possibly 
zero), outside this season. With the inclusion of data from 2002 and 2003, thus extending the 
Davis data base from 7 to 9 years, the winter effect changed to 4.8±1.3 K/lOOsfu [French et 
al, 2005], 

675 Offermann et al. [2004] have also reported different solar influence during different months 
of the year. The monthly responses suggest considerable variations even though the estimated 
error bars are large, because of the strong dynamical variability. The mean of the data is 3.0 ± 
1.6 K/lOOsfu. This is in agreement with the analysis based on annual means that gave 3.4 
K/lOOsfu. The authors concluded that long-term trend effects as measured at Wuppertal and 
680 solar cycle influences are almost statistically independent, which means that there is little 
interference between both types of results as noted by them. 

6. MODEL SIMULATIONS OF SOLAR RESPONSE 

685 

In principle, those models which are able to account properly for the vertical coupling 
processes in different altitudes of the atmosphere are suitable to study solar variability effects. 
There are only few model studies to assess the effect of solar variability on temperature or 
other parameters in the MLT region, in comparison to stratospheric regions where many 
690 models have been used [Rozanov et al., 2004 and references therein]. Some studies addressed 
the effect of the 27-day rotational variation with one-dimensional (1-D) [Brasseur et al., 
1987; Summers et al., 1990; Chen et al.. 1997] or 2-D [Zhu et al.. 2003] chemical dynamical 
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models. Current models of the effect of 11 -year solar cycle on the middle atmosphere 
temperatures are inconclusive. Studies by Brasseur [1993] and Matthes et al. [2004] suggest 
695 variations of upper mesospheric temperatures of about 2 K in response to changes in the solar 
UV flux, (red added extra in rev-1). The 1 1-yr solar cycle variability was studied with 
different versions of the SOCRATES (Simulation of Chemistry, Radiation, and Transport of 
Environmentally Important Species) interactive 2-D model by Huang and Brasseur [1993] 
and Khosravi et al. [2002], Huang and Brasseur [1993] arrived at a peak-to-peak temperature 
700 response to solar activity in the mesopause region of about 10K, whereas Khosravi et al. 
[2002] derived a value of 5K. Garcia et al. [1984] have reported a solar response of 6K 
between solar minimum and maximum activity using their 2-D model . 

Most of the initially developed general circulation models (GCMs) extend generally from the 
705 surface to the mid-stratosphere. Later, some of these GCMs have been extended to 
approximately 75-100 km altitude [e.g., Fels et al.. 1980; Boville. 1995; Hamilton et al., 
1995; Manzini et al., 1997; Beagley et al., 1997] or even up to the thermosphere [Miyahara et 
al., 1993; Fomichev et al., 2002; Sassi et al., 2002], Chemical transport models that treat 
chemical processes up to the mesosphere “offline” from the dynamics have also been 
710 developed [e.g., Chipperfield et al., 1993; Brasseur et al., 1997]. Coupled dynamical- 
chemical models covering this altitude range used mostly a mechanistic approach [e.g., Rose 
and Brasseur, 1989; Lefevre et al., 1994; Sonnemann et al., 1998] in which the complex 
processes of the troposphere are replaced by boundary conditions applied in the vicinity of 
the tropopause. However, all these three-dimensional upper atmospheric numerical models 
715 for the mesosphere and lower thermosphere usually do not include the troposphere. However, 
it is well known that mesospheric dynamics are largely determined by upward propagating 
waves of different kinds that have their origin, in general, in the troposphere. Only very 
recently, models have been developed which include a detailed dynamical description of the 
atmosphere including the troposphere, have their upper lid in the thermosphere, and can be 
720 coupled to comprehensive chemistry modules (GCMs with interactive chemistry are referred 
to as chemistry climate models, CCMs). Models of this type are the Extended Canadian 
Middle Atmosphere Model [EXCMAM, Fomichev et al., 2002], the Whole Atmosphere 
Community Climate Model [WACCM, e.g. Beres et al., 2005; Garcia et al., 2007], and the 
Hamburg Model of the Neutral and Ionized Atmosphere [HAMMONIA, Schmidt and 
725 Brasseur, 2006; Schmidt et al., 2006]. Only these recent models can be expected to 
realistically describe the atmospheric response to the variability of solar irradiance. 
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The newly developed HAMMONIA model combines the 3-D dynamics from the ECHAM5 
model with the MOZART-3 chemistry scheme and extends from Earth's surface up to about 
730 250 km. In the mesosphere and lower thermosphere the distance between the levels is 

constant in log-pressure and corresponds to about 2 to 3 km, depending on temperature. 
Schmidt et al. [2006] have performed model simulations on both the doubled CO 2 case and 
the role of the 11 -year solar cycle in trend studies. They find a temperature response to the 
solar cycle as 2 to 10 K in the mesopause region, with the largest value occurring slightly 
735 above the summer mesopause (ca. 100 km). Up to the mesopause, the temperature response 
may either be positive or negative depending upon longitude, particularly for middle and high 
latitude in winter. This study (like several other modeling reports) also points out the 
importance of distinguishing the presentation of results according to the choice of the vertical 
coordinate system, since the effects of subsidence look quite different at constant geometric 
740 than at constant pressure altitude. Marsh et al [2007] have recently used the WACCM 
(version 3) model. The response of the MLT region in WACCM3 is broadly similar to that of 
the HAMMONIA model shown in Schmidt et al. [2006]. Marsh et al. [2007] reported that the 
global-mean change in temperature for 50-80 km is between 0.3 and 1.5 K/lOOsfu, for 80-90 
km it ranges from 1.5 to 2.5 K/lOOsfu, and for 90-110 km it is 2.5 to about 
745 5 K/lOOsfu. There is a local minimum around 66 km above which the solar 

temperature response increases with increasing height. 

7. CONCLUSIONS 

750 As evident from the paper, the comparison of the results obtained by different observations 
separated by several decades is complicated. Nevertheless, there are a number of occasions 
where most of the temperature responses to solar variability indicate consistency and some of 
the differences are even understandable. The present status of MLT region solar response 
based on the available measurements can be broadly described as follows: 

755 

(1) Recognition of positive signal in the annual mean solar response of the MLT regions with 
an amplitude of a few degrees per lOOsfu. This agrees with numerical simulations of 
coupled models. 
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(2) Most Northern Hemispheric results indicate a solar response of the order of l-3K/100sfu 

760 near the OH airglow emission height at mid-latitude which becomes negligible near the 

Pole. 

(3) In the Southern Hemisphere, the few results reported so far indicate the existence of a 
stronger solar response near the pole and a weaker response at lower latitudes. 

(4) There is increasing evidence for a solar component of the order of 2-4K/100sfu in the 
765 lower thermosphere (92-100 km) which becomes negligible around 103 km. 

(5) In the mesosphere, the mid-latitude solar response of l-3K/100sfu is consistent with 
satellite, lidar and US-rocksonde data whereas Russian results indicate a more variable 
behavior. 

(6) Existence of a solar response of l-3K/100sfu for the mesospheric region in the tropics. 

770 

Most recent GCM model results indicate a higher upper limit of solar temperature response as 
compared to observations (2 to 10 K per solar cycle in the mesopause region), with the 
largest value occurrs slightly above the summer mesopause (ca. 100 km). Up to the 
mesopause, the temperature response may either be positive or negative depending upon 
775 longitude, particularly for middle and high latitude in winter. 

It is becoming increasingly evident that the solar response in the MLT region is highly 
seasonally dependent. This might explain the dispersion of the values in annually averaged 
solar response reported in this paper as it might have been caused by the seasonal distribution 
780 of observation. The differences in temperature response to solar activity in the mesopause 
region are mainly caused by changes in the vertical distribution of chemically active gases 
and by changes in UV irradiation. The intervention of dynamics (e.g., the mediation by 
planetary waves) further compounds the picture, which is only likely to become clearer, after 
more results on the long term solar response become available. Hence this topic remains a 
785 problem to explore more rigorously in the future. 

However, the major challenge is in the interpretation of the various reported results which are 
diverse and even indicate latitudinal and a considerable amount of longitudinal variability in 
solar response. The high degree of similarity in the response of the mesosphere to increasing 
790 surface concentration of greenhouse gases and to 11 -year solar flux variability suggests that 
climate change in the mesosphere may not be associated with anthropogenic perturbations, 
alone. If long-term increase in the well-mixed greenhouse gases, in particular C02, alters the 
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thermal structure and chemical composition of the mesosphere significantly and if these 
anthropogenic effects are of the same magnitude as the effects associated with the 11 -year 
795 solar cycle then the problem is more difficult to analyze. It is therefore necessary to 
discriminate between the two effects, and to identify their respective contribution to the 
thermal and chemical change in the mesosphere. The cyclic nature of the variability in solar 
UV flux over decadal time scales may provide the periodic signature in the observed response 
that could be used to identify variations in solar activity and other perturbations causing the 
800 changes, but this requires longer series of observations. 

8. OUTLOOK 

Looking forward, there are many compelling scientific questions and analyses still waiting to 
805 be addressed and undertaken. Amongst those relating to the long term response in the MLT 
region to solar variations are: 

• The detailed analysis of the trends of parameters like winds and minor constituent 
concentration (water vapor and ozone) in the mesopause region in order to properly 
understand the MLT temperature trends and solar response. 

810 • The monthly-to-seasonal long-term temperature trend and solar cycle response in the 

mesosphere, including the mesopause region. 

• Modeling studies of solar trends as derived from existing General Circulation Models 
(GCMs). The consistency between observed and modeled temperature, radiation, and 
chemistry must be evaluated. These studies are expected to yield future measurement 

815 recommendations. 

Finally, we expect that great progress in understanding the MLT response to solar variations 
will be provided by the NASA TIMED mission that has just completed 6 years in orbit. 
Strong evidence for solar cycle influence on the infrared cooling of the thermosphere has 
820 already been shown by Mlynczak et al., [2007]. They noted a factor of 3 decrease in the 
power radiated by NO in the thermosphere from the start of the mission (near solar 
maximum) through calendar year 2006, which corresponds nearly to solar minimum. The 
TIMED data set with its measurements of temperatures, constituents (ozone, water vapor, 
carbon dioxide, O/N 2 ratio, etc.) and solar irradiance will enable a unique dataset from which 
825 the effects of the 11 -year solar cycle can be confidently determined. Efforts are now 
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underway to secure operation of the TIMED mission through the next solar maximum in 
approximately 4-5 years. 
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FIGURE CAPTIONS: 


1170 

Figure 1: Ground based mesospheric and lower thermospheric data stations for the 
measurements of temperature. The FPI represents Fabry-Perot Interferometer and ISR 
represents Incoherent Scatter Radar. Those stations whose results have been used in the 

present paper are marked with *sign. 

1175 

Figure 2: Atmospheric temperature response to the 27-days and 11 -years solar cycle as a 
function of altitude from different studies focus on tropical regions. SAMS [Keating et al., 
1987] and PMR [Ebel et al., 1986] analyses are performed on the solar rotation time scale. 
HALOE analyses above the ±20° latitude band [Remsberg, 2007] concern the 1 1 -years time 
1180 scale. Two different analyses [Beig and Fadvanis, 2003, Kokin et al., 1990] have been used 
to analyze the soviet rocket station of Thumba. 

Figure 3: Atmospheric temperature response to the 11 -years solar cycle as a function of 
altitude from different studies covering mid-latitudes. Error bars associated with the OHP 
1185 lidar correspond to the range of seasonal changes [Keckhut et al., 2005]. The hatched area 
corresponds to the HALOE observations in different latitude bands [Remsberg, 2007]. ‘US 
Rocket’ correspond to an average of 3 subtropical rocket stations [Keckhut et al., 1999]. The 
Soviet rocket sites (Volgograd, Molodezhnaya, Heiss) have been analysed by Kokin et al. 
[1990] and Golitsyn etal. (2006). 

1190 

Figure 4: Solar cycle response in the Northern Hemisphere for the mesopause region. 

Figure 5: Solar cycle response in the Southern Hemisphere for the mesopause region. 

Figure 6: Solar cycle response in the lower thermospheric region. 

1195 Table 1. Solar heating, chemical heating, and radiative cooling terms, and associated airglow 
losses, for the mesosphere and lower thermosphere. 

Table 2. Solar response in Temperature in the Mesosphere (50-79 km) Region 
1200 Table 3. Solar response in temperature of the mesopause region 
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1205 


1210 


1215 


1220 


Table 1. Solar heating, chemical heating, and radiative cooling terms, and associated airglow losses, 
for the mesosphere and lower thermosphere 

Solar Heating 

Chemical Heating 

Radiative Cooling 

0 2 Schumann- Run ge Continuum 

H + 0 3 -» OH + 0 2 

C0 2 (v 2 ) - 15 nm 

0 2 Schumann-Runge Bands 

H + 0 2 + M H0 2 + M 

0 3 (v 3 ) - 9.6 nm 

0 2 Lyman-alpha 

0 + OH H + 0 2 

NO(o) - 5.3 nm 

0 2 Atmospheric Bands 

o+ho 2 -»oh + o 2 

H 2 0 rotational - 25-50 nm 

0 3 Hartley Band 

o + o + m-»o 2 + m 

0( 3 P) - 63 nm 

All radiation X < 120 nm 

0 + 0 2 + M O 3 + M 



0 + O 3 + M 0 2 + 0 2 


Associated airzlow losses 



0 2 ('A) - 1.27 nm 

OH(u = 1 to 9) 


0 2 ('E) - 762 nm 

OzC'A) - 1 .27 nm 


C0 2 (v 3 ) - 4.3 pm 

03 ^ 3 ) - 9.6 pm 
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Table 2. Solar response in temperature in the mesosphere (50-79 km) Region 

1225 


Reference 

Technique 

Years of 
Analysis 

Location 

Height (km) 

11 year Solar cycle 
response (K) 

Remarks 

Keating et al. 

SAMS 

1980’s 

tropics 

70 

2.5 

27 days 

(1987) 




60 

1.2 






50 

1.1 


Ebel et al., (1986) 

PMR 

1975-1978 

tropics 

50 

3.5 

27 days 





60 

3 






70 

2.3 


Remsberg, (2007) 

HALOE 

1991-2005 

tropics 

50 

0.9 






60 

1 






70 

1.2 


Remsberg, (2007) 

HALOE 

1991-2005 

mid-latitude 

50 

0.95 






60 

1.5 






70 

2.05 


Beig and Fadnavis 

Indian- 

1971-1993 

8.5°N, 77°E 

50 

-0.8 


(2002) 

rocketsonde 



60 

0.2 






70 

3 


Kokin et al., (1990) 

Indian- 


8,5°N, 77°E 

50 

-0.5 



rocketsonde 

1971-1993 


60 

1.8 






70 

1.2 


Keckhut et al., (2005) 

us- 

1969-1995 

22°N-34°N 

50 

1.8 



rocketsonde 



60 

2.4 






70 

0.4 


Keckhut et al., (2005) 

Lidar 

1979-1997 

44°N, 6°E 

50 

0.75 






60 

2.2 






70 

1.7 


Golitsyn et al. (2006) 

Volgograd 

1971-1993 

48.7°N, 44°E 

50 

-1.8 


and Kokin et al. 

rocketsonde 



60 

4 


(1990) 




70 

3 


Golitsyn et al. (2006) 

Heiss Island 

1971-1993 

80.6°N, 58°E 

50 

-3 


and Kokin et al.(1990) 

rocketsonde 



60 

3 






70 

5 
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Golitsyn et al. (2006) 
and Kokin et al.(1990) 

Molodezhnaya 

rocketsonde 

1971-1993 

68° S, 46°E 

50 

60 

70 

1 

0.5 

5 


Liibken (2000a, b) 

Rocket 
grenades 
Falling spheres 

1960s, 

1970+ 

1987-2000 

69°N, 10°E 

50-80 

no response 
(summer) 

most observations 
during solar minimum 
1960s, 80s, 90s 
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Table 3. Solar response in temperature of the mesopause region 


(a) Lower mesopause region (80-90 km) 

Reference 

Technique 

Years of 
Analysis 

Location 

Height (km) 

IjjlEMMB 

Remarks 

NORTHERN 

HEMISPHERE 







Sigemes et al. 
(2003) 

Hydroxyl rotational 
temperature 
Bands: 6-2 
Bands: 8-3 

1983-2002 

1980-1983 

78°N, 15°E 

OH emission 
layer 

+0.0 (±0.5) (winter) 


Liibken (2000, 
2001) 

Rocket grenades 
Falling spheres 

1960s, 

1970+ 

1987-2000 

69°N, 10°E 

80 

no response (summer) 

most observations during 
solar minimum 1960s, 80s, 
90s 

Espy and Stegman 
(2002) 

Hydroxyl rotational 
temperature 
Band: 3-1 
also 4-2 

1991-1998 

59.5°N, 18°E 

OH emission 
layer 

1.6 (±0.8) (winter) 

consistent with +2 K/lOOsfu 
[Espy, personal 
communication 2007] 

Pertsev and 
Perminov (2007) 

Hydroxyl rotational 
temperature 
Band: 6-2 

2000-2006 

55.7°N, 37°E 

OH emission 
layer 

7.5 (±0.5) (winter) 
2 (±2.0) (summer) 

4.5 (±0.5) (annual) 

Also intensity response for 
OH and 0 2 bands 

Offermann et al. 
(2004) 

Hydroxyl rotational 
temperature 
Band: 3-1 

1980-1984 

1987-2002 

51.3°N, 7°E 

99 

+3.0 (±1.6) 


Lowe (2002) 

Hydroxyl rotational 
temperature 
Band: 3-1, 4-2 

1989-2001 

42.5°N, 81°W 

99 

+1.5 (± 1.1) 


Golitsyn et al. 
(2006) 

Combined (Major: 
Hydroxyl rotational 
temperature 
Bands: several) 

1976-1991 

41.8°N, 43° E 
and 

55.7°N, 37°E 

80-95 

-5 (±1.7) (winter) 
+8 (±1.7) (summer) 
+2 .7 (± 1 .7) (annual) 


She and Krueger 

Na-Lidar 

1990-2001 

41 °N, 105°W 

83 


update C.-Y.She, personal 
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(2004) 

Remsberg, and 
Deaver (2005) 


UARS -HALOE - 
sunrise, sunset 
(satellite) 


1991-2004 


40° S - 40°N 


87 

80 


SOUTHERN 

HEMISPHERE 


Azeem et al. 
(2007) 

Hydroxyl rotational 
temperature 
Bands: 3-1 

1994-2004 

90° S 

Antarctica 

OH-emission 
layer height (87) 

Hernandez et al. 
(2003) 

Hydroxyl rotational 

temperature 

Bands:6-2 

1991-2002 

90° S 

99 

French et 
al.(2005) 

Hydroxyl rotational 
temperature, 

Band: 6-2 

1990, 

1995-2003 

69° S, 78°E 

99 

French and Bums 
(2004) 

Hydroxyl rotational 
temperature, 

Band: 6-2 

1995-2001 

69° S, 78°E 

99 

Reisin and Scheer 
(2002) 

Hydroxyl rotational 
temperature 
Band: 6-2 

1986, 1987, 
1992, 1997- 
2001 

32° S, 69° W 

99 

Scheer et al. 
(2005) 

Hydroxyl rotational 
temperature 
Band: 6-2 

1998-2002 

32° S, 69° W 

99 

Clemesha et al. 
(2005) 

Hydroxyl rotational 
temperature 
Band: 9-4 
Bands: 6-2 

1987-1994 

1994-2000 

23°S, 45°W 

99 


(b) upper mesopause region and lower thermosphere (90-110 km) 


Reference 

Technique 

Years of 
Analysis 

Location 

Height (km) 

Golitsyn et al. 

estimated from 

1976-1991 

41.8°N, 43° E 

92 


3 (±1) 

communication [2007] 

0.5-1. 7 amplitudes 

range of solar response 

(±0.6)-> 1.5 (±.0.8) 

estimated from published 

K/lOOsfu 

amplitude 



4.0 (±1.0) 

(winters) 

13.2 (+3.5) (winter) Doppler (kinetic) 

temperature 

4.8 (±1.3) (winter) 

6-7 (mid-winter) published range 

no response 


0.92 (±0.32) 


8-10 


range of solar response 
estimated from published 
amplitude 
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(2006) 

Sodium emission 
intensity 


and 

55.7°N, 37°E 




Golitsyn et al. 
(2006) 

estimated from 
Atomic oxygen 
emission intensity 

1976-1991 

41.8°N, 43° E 
and 

55.7°N, 37°E 

97 

2.3 (+1.1) 


She and Krueger 
(2004) 

Na-Lidar 

1990-2001 

41 °N, 105°W 

92-98 

103 

4.0 (±1.8) 

0.0 (±1.4) 


Reisin and Scheer 
(2002) 

0 2 (0-1) rotational 
temperature 

1986, 1987, 
1992, 1997- 
2001 

32° S, 69° W 

0 2 emission 
height (95 km) 

1. 6-5.6 

range of upper limit 
estimates 

Scheer et al. 
(2005) 

0 2 (0-1) rotational 
temperature 

1998-2002 

32° S, 69° W 

O 2 emission 
height (95km) 

3.3 (± 0.3) no trend 1.32 
(± 0.3) with temporal trend 

mean response. See text for 
details 
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Figure 1 : Ground based mesospheric and lower thermospheric data stations for the measurements of temperature. The FPI represents Fabry-Perot 
Interferometer and ISR represents Incoherent Scatter Radar. Those stations whose results have been used in the present paper are marked with *sign. 



1235 



Altitude -km 

- - ■ SAMS-27d PMR-27d 

— - Thumba-lly-Beig HALOE-lly 
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Figure 2: Atmospheric temperature response to the 27-days and 11 -years solar cycle as a 
function of altitude from different studies focus on tropical regions. SAMS [Keating et al., 
1987] and PMR [Ebel et al., 1986] analyses are performed on the solar rotation time scale. 
HALOE analyses above the ±20° latitude band [Remsberg, 2007] concern the 11 -years time 
scale. Two different analyses [Beig and Fadvanis, 2003, Kokin et al., 1990] have been used to 
analyze the soviet rocket station of Thumba. 
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Figure 3: Atmospheric temperature response to the 1 1 -years solar cycle as a function of altitude 
from different studies covering mid-latitudes. Error bars associated with the OHP lidar 
correspond to the range of seasonal changes [Keckhut et al., 2005]. The hatched area 
corresponds to the HALOE observations in different latitude bands [Remsberg, 2007]. ‘US 
Rocket’ correspond to an average of 3 subtropical rocket stations [Keckhut et al., 1999]. The 
Soviet rocket sites (Volgograd, Molodezhnaya, Heiss) have been analysed by Kokin et al. [1990] 
and Golitsyn et al. (2006). 
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Figure 4: Solar cycle response in the Northern Hemisphere for the mesopause region. 
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Figure 5: Solar cycle response in the Southern Hemisphere for the mesopause region. 
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Figure 6: Solar cycle response in the lower thermospheric region. 
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